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The evolution of the Kolmogorov-Zakharov (K-Z) spectrum of weak turbulence is studied in the 
limit of strongly local interactions where the usual kinetic equation, describing the time evolution 
of the spectral wave-action density, can be approximated by a PDE. If the wave action is initially 
compactly supported in frequency space, it is then redistributed by resonant interactions producing 
the usual direct and inverse cascades, leading to the formation of the K-Z spectra. The emphasis 
here is on the direct cascade. The evolution proceeds by the formation of a self-similar front which 
propagates to the right leaving a quasi-stationary state in its wake. This front is sharp in the sense 
that the solution remains compactly supported until it reaches infinity. If the energy spectrum has 
infinite capacity, the front takes infinite time to reach infinite frequency and leaves the K-Z spectrum 
in its wake. On the other hand, if the energy spectrum has finite capacity, the front reaches infinity 
within a finite time, t* , and the wake is steeper than the K-Z spectrum. For this case, the K-Z 
spectrum is set up from the right after the front reaches infinity. The slope of the solution in the 
wake can be related to the speed of propagation of the front. It is shown that the anomalous slope 
in the finite capacity case corresponds to the unique front speed which ensures that the front tip 
contains a finite amount of energy as the connection to infinity is made. We also introduce, for 
the first time, the notion of entropy production in wave turbulence and show how it evolves as the 
system approaches the stationary K-Z spectrum. 

PACS numbers: 04.30.Nk, 47.35.-|-i, 92.10.Hm 



I. INTRODUCTION AND MOTIVATION 



Wave turbulence is concerned with the statistical description of an infinite sea of dispersive waves, which are 
weakly coupled by nonlinear interactions and maintained away from equilibrium by interaction with sources and sinks 
of energy. The theory has found practical application in many branches of physics including the description of surface 
waves on fluid interfaces H 011 m, Alfven wave turbulence in astrophysical plasmas nonlinear optics 3 and 

acoustics |lO|,|l^ to name a few. 

The central quantity of theoretical interest is the spectral wave action density, nk, which describes how the exci- 
tations in the system are distributed among different wave- vectors, k. Under fairly weak assumptions |ll|, the long 
time behaviour of rik is given by an equation known as the wave kinetic equation. For a system dominated by four 
wave interactions this equation takes the form, 



dt 



47r J iTkkik.ksl ' .f4K]5(k+ki-k2-k3) dkidkadkg, (1) 



where 



[uk] = nkHkiTik^nka ( — + — —] (5(wk+t^ki-t^k2-'^k3)- (2) 

\nk rtki "ka "-k:,/ 

Equation Q is the analogue for waves of the Boltzmann equation of classical statistical mechanics. In many applica- 
tions, Wk and Tkkikaka are homogeneous functions of their arguments. Their degrees of homogeneity shall be denoted 
by a and 7 respectively. Under rescaling, k Ak, they transform as follows : 

cjAk = A"tJk (3) 

^AkAkiAkaAka ^'''Tkkikakg • (4) 

It was shown by ZakharovTsj in the 60's that if the energy sources and sinks are separated by an "inertial range", 
equation |^ has exact isotropic steady state solutions, 

nk = CiPl/3fc-(27+3rf)/3 (5) 

rik = C2Q'/'k-^^^+"'-''>/\ (6) 
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which carry constant fluxes of conserved densities, in this case energy flux, P, or wave action flux, Q, between sources 
and sinks. These steady state spectra are the direct analogues of the direct and inverse cascades in hydrodynamic 
turbulence and are referred to as Kolmogorov-Zakharov (K-Z) spectra. They have been well observed experimentally 
in a variety of contexts. 

In 1991 Falkovich and Shafarenko 5] addressed the question of how the K-Z spectrum is set up in time if rtk is 
initially compactly supported in wave-vector space. They used a self-similar solution of |^ and an assumption that, 
for the direct cascade, the total energy increases linearly in time, to show that (|3J) is set up by a nonlinear front which 
propagates towards k = oo and leaves the /j-(27-i-3rf)/3 spectrum in its wake. 

However subsequent numerical simulations of the kinetic equation for Alfven wave turbulence performed by Galtier 
et al. 1^ suggested that the development of the K-Z spectrum may proceed by a different route. For the Alfven wave 
system, the K-Z energy spectrum has finite energy capacity, meaning that on the K-Z spectrum, Ji?(k) dk < oo. 
This implies that the nonlinear front must reach A: = oo within a finite time, t*. They noticed that the spectrum in 
the wake of the front was significantly steeper than the K-Z value for times less than the singular time, t* , and that 
the K-Z spectrum then developed from right to left after the front reached k ^ oo. Other work by Pomeau et al.0 on 
the inverse cascade in the Nonlinear Schrodinger equation suggested that there might be anomalous quasi-stationary 
spectra associated with non-stationary solutions of kinetic equations. However no-one has yet made a specific attempt 
to search for them. 

One of the challenges in far-from-cquilibrium systems is to understand the means by which stationary states are 
reached and to ask if there are functionals analogous to the entropy in equilibrium systems. What we will show is 
that while the entropy, which for wave turbulence is formally (see for example 0) 



is not well defined on the steady state solutions, spectra (|3J) and 10, its production rate is. We find that for < t < t*, 
when the spectrum in the wake of the front is steeper than the K-Z spectrum, the entropy production is positive. At 
t* the connection to k = oo is made and energy is no longer a conserved quantity. For t > t*, the K-Z spectrum is 
established via a front which travels back from k = oo. During this stage, the entropy production rate, while still 
positive, gradually decreases and asymptotes to zero, its value on the exact K-Z spectrum. We conjecture that this 
scenario, established in this paper for the differential approximation to the kinetic equation, will be widely valid 
for finite capacity non-equilibrium systems, including three-dimensional hydrodynamic turbulence at large Reynolds 
numbers. 

This leads us to the topic of this article. We have made an extensive study of the non-stationary solutions of the 
so-called differential kinetic equation of local wave turbulence. This model equation is obtained from ^ under the 
assumption that the interaction co-efficient, Ikkik2k3 is strongly localised in kkik2k3 space. It has the advantage of 
replacing the integro-differential kinetic equation with a PDE. We find that the qualitative bahaviour observed by 
Galtier et al is present in this model in the finite capacity case. Since we are dealing with a PDE, we can go a lot 
further in terms of understanding. 

The organisation of the article is as follows. In section^lwe introduce the differential kinetic equation and describe 
a few of its properties which make it a good model of wave turbulence. We also introduce exact expressions for 
the fluxes of energy (P), the flux of particles (Q), and the entropy production rate in terms of its flux (i?) and 
its bulk production rate (T). The latter is always positive definite. We calculate each of these quantities on the 
algebraic solutions, oc fc~"^. Section HTll contains the details of some numerical simulations of the PDE. These 
simulations suggest that the nonlinear front is "sharp" in the sense that remains compactly supported for t < t*. 
There is a singularity, a divergence in the second derivative in fact, at the front tip between the regions nk = and 
Uk > 0. Next, in section Hvl we construct a family of self-similar solutions of the differential kinetic equation which 
are parameterised by a single free parameter. This free parameter can be interpreted as the asymptotic slope behind 
the front. Following that, in sectional we use this self-similarity analysis to formulate a hypothesis which we call the 
critical front speed hypothesis. This hypothesis is based on physical arguments and allows us to select a critical value, 
Xc, for the asymptotic slope, given by 



where Xq denotes the usual K-Z exponent for the direct cascade. This formula is well supported by our numerical 
simulations. In our conclusion, we attempt to make a connection with entropy production arguments. Two appendices 
are provided. In the first, we analyse the mathematical structure of the similarity equation and try to understand how 
the critical slope is related to the solution trajectories of the underlying ordinary differential equation. The second 
appendix contains a brief outline of the numerical methods used. 




(7) 



Xc = Xo + 



27- 3a 
12a 



(8) 
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II. THE DIFFERENTIAL KINETIC EQUATION 



We begin by briefly discussing the origin of the differential kinetic equation. Assuming that the wave-action 
spectrum rapidly becomes isotropic, and averaging over angles, we can make a transformation from d-dimensional 
wave- vector space to frequency space, 



dt 

where 



Suuii^2i^3 J^A [nj\ dujiduj2duj3 (9) 



= ^t: f\Tkk,k,k,\ '(5(k+ki-k2-k3) (kk.hhf-'^^^^ dn (10) 
J aoj dull aoj2 aojs 

and is defined by requiring that 

' (t>{iu)N^duj = I <^(|kr) nkdk, (11) 



for any test function, cj). Here the volume element dU, represents integration over the angular variables in kik2k3 
space and the wave-vector moduli are related to the frequency via the dispersion relation, 

Uk = cfc«. (12) 

If we assume that the interaction coefficient, Tkkikaka is strongly local in kkik2k3 space then can be approximated 
by a differential equation 

dN^ . 92 ^ . . / 1 

where 



n^{t) = n{k{uj),t) (14) 
s = 3x0 -f 2 (15) 

= — 5 ■ (16) 

/ is a pure number which comes from the angular integrations in (|10|) . This equation is called the differential kinetic 
equation. 

The local approximation leading to the differential kinetic equation, \Vi\ . is rather drastic and is not justified for 
many of the physical applications of weak turbulence. Nonetheless, it retains many of the qualitative features of the 
full kinetic equation, jnj. It provides an excellent model in the context of which these features can be easily studied. 
In particular, the differential kinetic equation respects the conservation laws embodied within its integro-differential 
precursor, namely the conservation of the total energy. 



E = / ljN^ duj, (17) 

and the total number of particles, 

N = I duj. (18) 



It also preserves the scaling and homogeneity properties of the kinetic equation. Consequently, the pure scaling 
solutions of the kinetic equation, the equilibrium thermodynamic spectra and the non-equilibrium K-Z spectra, are 
also solutions of the differential kinetic equation. These solutions are of the form = clu^^ where x takes one of the 
following values: 

a; = 1, a; = 0, (thermodynamic) 

or 

_ 27 + 3rf _ 2j + 3d-a 

3^ ' ^^-^>- 
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FIG. 1: Plot of the bulk entropy production, T, the entropy flux, R, and the total entropy production, S, as a function of the 
exponent, x, of the spectrum for power law spectra, riu} — cw~^ . 



It is convenient to define 



The two conservation laws can be written as continuity equations : 

dN^ dQ 



IC[n^]=u;'^o+'nt-^(—]. (19) 



dt dijj 
dE^ dP 



dt duj 



= (20) 
= 0, (21) 



where 



is the local flux of particles, 



is the local flux of energy and 



P = JC~Lo^ (23) 

OUJ 



^ LuN^ (24) 



is the energy density if frequency space. Note that P is defined to be positive when energy flows to the right in lo 
space and Q is positive when particles flow to the left. 

The entropy, S, of a wave system is formally J Inrtk dk. The production rate. 



can readily be calculated to be 



- = - /lnn,dk= /-— d.., (25) 



S^-^+T. (26) 
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The entropy flux, i?, which is positive for entropy flow to large wave- numbers, is 



Q d P dn^ 
5--5-(w?^a;) 5--H— (27) 



and the local bulk entropy production rate, T, is 



T = /..-o+2„4 ( f -1 ) ) . (28) 



Note that T is positive definite and zero on the thermodynamic solution, = t / (uj — /j.) , r being the temperature 
and /i the chemical potential. Indeed if the system were isolated, say on the interval lui < lu < uj2, and P, Q and R 
were identically zero, then the usual principles of equilibrium systems would apply. 

On the solution = cuj~^, the quantities K, Q, P, R, T and — ^ + T are calculated. They are: 

K = Ic^uj^''°~^''x{x-l) (29) 

Q ^ 3Ic^x(x~1)(xq-x)uj^''°-^'=-^ (30) 

P = 3Ic^x{x-l){x- {xo- ^))uj^''°-^'' (31) 
T ^ c2/a;2(x- 1)2^^3.0-2.-2 (32) 

R = 3c2/x(a;- l)(xo - |a;)a;3^''-2^~\ (33) 
^ = -^+r = 9c2/a;(x-l)(x-a;o)(a;-(xo-i))w3^"~'"-^ (34) 

Similar expressions, which have the same zeros (as functions of x) obtain for the case when the differential approx- 
imation is replaced by the full collision integral. In particular, we note that the entropy production rate is always 
positive (we assume Xq > 4/3) for 1 < x < xq — 1/3 and for x > xg, the K-Z exponents for particles and energy 
respectively. The relevant functions are plotted in figure ^ For xq — 1/3 < x < xq, the entropy production rate is 
negative. This corresponds to a situation when the particle flux is building a condensate state |9(|. 



III. NUMERICAL OBSERVATIONS OF NON-STATIONARY SPECTRA 



Let us now turn our attention to non-stationary solutions of (|13|l . In particular we are interested in how the K-Z 
spectra are set up if we begin from an initial condition which is compactly supported at low frequencies. Early work 
on this question focussing on the direct cascade by Falkovich and Shafarenko |5| suggested that the K-Z spectrum is 
set up by a nonlinear front which propagates to the right leaving the K-Z spectrum in its wake. Recent numerical 
studies by Galtier et al suggest that this problem is more subtle. They found that in the case where the K-Z 
energy spectrum has finite capacity, the approach to the steady state proceeds by a different mechanism. For this 
system the nonlinear front reaches infinite frequency within a finite time, t*. They found that the quasi-stationary 
spectrum in the wake of this front was actually steeper than the K-Z spectrum. The K-Z spectrum was then set up 
from right to left after the front reached infinity. 

We investigated whether there was evidence of this behaviour in the solutions of the differential kinetic equation. 
We solved the differential kinetic equation numerically and followed the evolution from an initial condition compactly 
supported at low frequencies. The results presented in this section were obtained by allowing the initial data to decay 
freely in the absence of forcing or damping. Some details of the numerical methods used are contained in appendix 

m 

Let us first consider the situation in which the energy spectrum has infinite capacity, 27 < 3a. Figure 12 shows a 
sequence of snapshots of and the local slope at successive times for a system with parameter values a = 0.5,7 = 
0, d = 1. We conclude that the slope tends to the K-Z value far behind the front in agreement with previous 
expectations. The K-Z spectrum is set up asymptotically in time from the left of the window of transparency to the 
right. 

Now consider what happens when we increase the value of 7 and bring the system into the finite capacity regime. 
Figure 13 shows the absolute value of the local slope of at similar stages in the evolution for a sequence of values 
of 7 in the finite capacity regime. 7 takes the values 2.0, 2.5, 3.0 and 3.5. Visually, it is clear that as 7 increases, 
the slope of the spectrum behind the front tends to a value which is increasingly steeper than the K-Z value. This 
observation is supported by fitting power law functions to the numerical data. The "best fit" slopes are presented in 
Table ^ along with the K-Z values for comparison. The slope given by equation © agrees well with the numerical 
observations. The time, t*, required for the front to reach infinity depends both on the parameters 7 and a and on 
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FIG. 2: Time evolution of the spectrum and the local slope for the parameter values a — 0.5,7 = 0,d = 1. The logarithmic 
scales are to base 10. 




FIG. 3; Snapshots of the non-stationary spectrum for a range of different values of the parameter 7, keeping a — 1/2 and 
d — 1. The logarithmic scales are to base 10. 
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xkz 


Xnum 




2.0 


-4.67 


-5.12 


-5.08 


2.5 


-5.33 


-5.93 


-5.92 


3.0 


-6.00 


-6.67 


-6.75 


3.5 


-6.67 


-7.50 


-7.58 



TABLE I: Numerical evidence for the anomaly. The table shows the values of the wake slopes obtained by fitting numerical 
data for a range of values of 7 in the finite capacity regime. These are to be compared to the K-Z values and the values of the 
critical slope obtained from equation Q : Xc = xkz + (27 — 3a)/12Q. 
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FIG. 4: Log-log plot showing the initial stages of relaxation to the K-Z slope after the front reaches uia. Parameter values are 
7 = 7/2, a = 1/2 and d = 1. The absolute value of the local slope is plotted as a function of frequency for three different times, 
just before t* , approximately at t* and after t* . There is a clear transition from the steeper slope, Xc, to the shallower xkz- 



the initial energy distribution. We have not, in the present work, made any systematic attempt to understand this 
aspect of the problem. 

The next obvious question is that of how the system makes the transition from this quasi-steady anomalous regime 
to the K-Z spectrum which we know to be the final steady state. For finite capacity systems this transition begins 
once the front reaches infinity. Clearly we cannot easily treat the divergence of the front numerically. Instead we allow 
the front to propagate into a regime of strong damping to the right of the window of transparency. This damping 
region is intended to mimic the energy sink role provided by the point uj — go. Figure^ shows the local slope of 
for successive times after the front reaches the dissipation region. The results of Figure^lare for 7 = 3.5 where there 
is a significant difference of about 0.9 between the anomalous slope and the K-Z slope. We see that the system begins 
to relax towards the K-Z spectrum as soon as the dissipation scale is reached. Notice that the relaxation is occurring 
from right to left. 

Another issue which is crucial to our explanation of the finite capacity anomaly is the structure of the front itself. 
From our numerical simulations we found that the front tip seems to be "sharp" in the sense that the spectral wave 
action density remains on compact support during the time evolution. Such sharp fronts have been known to exist 
in solutions of certain classes of nonlinear diffusion equations dating back to the work of Zel'dovich in the 50's. 
Subsequent work by Lacey et al. 0| showed that such fronts can be stationary, moving or exhibit waiting time 
behaviour where the tip remains stationary for a finite time before beginning to move. 

Let us suppose that our system possesses a sharp tip located at and 



n{LO, t) 



= for cj > uj*{t) 
> for cj < uj*{t). 



To the left of we look for a solution of equation H13() in the form of a power series, 

n{LO,t)^ A{t)^a,{Lu*{t)-u)'^\ (35) 

i 

where the constants ai and the exponents are to be determined by expanding equation {TJJ around uj — ^*(^) and 
substituting this representation for n{uj^t). The leading order term on the LHS is 

Aw*°"Voao(w* -w)™°-i-^. (36) 

The leading order term on the RHS is 

A^Lu*'mo{mQ + 1)(3too - 2)(3mo - 3)ag(w* - ojf'"-"-^. (37) 
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FIG. 5: Inset shows a detailed view of the structure of the front tip seen in the evolution of a finite capacity spectrum. The 
parameter values are 7 = 7/2, a = 1/2 and d — 1. The line is a least squares fit of the numerical data to the functional form 
n{u!) ~ A{u>* — (jj)'^ as suggested by equation 13911 . The fitted value of c is 1.52. 

Comparing powers of u* — lu immediately yields mo ~ 3/2. We cannot fix the time dependence although wc get the 
following relation between A{t) and uj*{t), 

Lo*{t)i-'^ = A\t)u.*{tY, (38) 

after choosing Oq ^ — {rtiQ + l)(3mo — 2)(3too — 3). The presence of a sharp tip and the structure 

n{uj,t) ^ A{t)ao{uj*{t)^ujf/'^, (39) 
immediately behind the tip is well supported by our numerical simulations as shown in figure |S1 

IV. SELF-SIMILAR SOLUTIONS OF THE DIFFERENTIAL KINETIC EQUATION 

To study non-stationary solutions of H13|) analytically, make the following self-similarity ansatz for the form of the 
solution, 

n^{t)^u:*{t)-- F(i^), (40) 

where the self-similar variable, rj, is defined by 

^ = ^y (41) 
Under this change of variables the derivatives transform according to the relations. 



dt \duj* di] 



and equation (|13|) can be rewritten as 



^,i-.,-.F-„i£).Jl(„-^.|,(l,,, (44, 
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where 

A = c^*'"+''/"-^"''. (45) 

Here xq = (27 + 3d)/3a is the exponent of the K-Z energy spectrum. The free parameter, x, is the asymptotic slope 
of the spectrum far behind the front. This is because we assume that there exists a quasi-stationary regime far behind 
the front, which is a simple power law to leading order, ~ uj"'^ . We must therefore have y — x in order to cancel 
the time dependence from the leading order part of (|40() . 

For a self-similar solution, A must be time independent. We are interested in situations where the system has finite 
energy capacity and generates a singularity within a finite time which we shall denote by t* . The appropriate solution 
of H45|l describing such situations is 

uj*{t)^{t* -t)K (46) 

It is convenient to define 

= — 5 ■ (47) 

The direct cascade has infinite energy capacity for < and finite energy capacity for Kd > 0. Upon substitution of 
the form ll^ into it follows that 

b =. {2ix ~ xo) - Ka)-' (48) 
A = -b (49) 

We can also define 

27 — a 



3a 



(50) 



such that the inverse cascade has infinite particle capacity for Ki > and finite particle capacity for Ki < 0. Equation 
(|48|l can be written in the equivalent form 

b^ (2{x-xo + 1/3) ~ K^y^ (51) 

which is more appropriate for studying the inverse cascade. Notice that the speed of propagation of the front at to*, 
as measured by b, is related to the asymptotic slope, x, behind the front. 

If we are interested in the direct cascade then — > 00 as i t*. This corresponds to & < 0. Conversely, for the 
inverse cascade, uj* ^ a.s t ^ t* , corresponding to 6 > 0. 

Upon substitution of (gHIl and g^l) into we obtain 

,2(. - _ {.F + ,f ) = ^ (1)) - (52) 

Let us now consider the structure of the front tip in the self-similar variables. From our numerical simulations we 
expect that there is a singularity in the solution as 77 ^ 1. For the direct cascade, we look for a singularity of the form 
(1 — 77)™, approaching from the left. For the inverse cascade we expect a singularity of the form (77— 1)™ approaching 
from the right. Let us restrict our attention to the direct cascade. 

Consider an expansion of the solution to the left of 77 = 1 in the form 

00 

F(77) = (l-77)"^a„(l-77)" (53) 

Tl = 

Taylor expand equation (|44|1 to the left of 77 = 1 and substitute this expansion for F. Matching of the leading order 
divergences fixes 777, = 3/2. In principle the coefficients, a„, can be computed to arbitrary order by matching powers 
of 1 — 77. The first few of them, computed using Mathematica, are 

ao = ^^3^/^ (54) 

ai = {—3d + 7as~2ax) (55) 

27 v6 a 

02 = — (18405 d^-Garf (4860 + 7931s - 3490 a;) 

1075032 V6a2 ^ 

-ha^ (47117 5^- 20 a; (1944 359 a;) - 28 s (-4374 + 833 a;))) (56) 
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FIG. 6: Application of the similarity transformation, (1401 . to the direct cascade at three consecutive times We take a range 
of values for the parameter x. Graph A takes x = 7.66, graph B takes x = 7.48, graph C takes x = 7.25. The self-similarity 
of the direct cascade and the front tip is evident from graph B. For reference, graph D takes x = 6.66, corresponding to K-Z 
scaling, and uses a log plot to display the entire function F(rj). It is clearly not self-similar for this value of x. 



So far, the slope behind the front, x, or equivalently by H48|) . the front speed, 5, is a free parameter in the similarity 
transformation. We wish to understand how the system picks the particular value for this parameter. If we are to 
observe the Kolmogorov-Zakharov slope behind the front then b must have the value 

bKZ = -l/^d, (57) 

which follows from putting x = xq in equation H48() . In section IIIII we presented numerical evidence that the slope 
behind the front was anomalous. We can also confirm this anomaly and check the correctness of the self-similarity 
argument outlined in this section by applying the transformations 14()|l . H41|l to the solution of the PDE. In figureElwe 
have taken a particular solution of l|13|) at a number of successive times and applied the similarity transformation for 
a selection of values of the free parameter x. The data presented is for the case a = 1/2, d = 1 and 7 = 7/2. It is clear 
that the direct cascade is self-similar for x — 7.48 approximately. This is to be compared with the theoretical prediction 
Xc — 7.56. The corresponding transformation with K-Z value, x — 6.66 in this case, is clearly not self-similar. 

The behaviour of the solution of the equation for the spectrum after the singular time, t* , is an interesting question. 
Let us consider what could be called the "bouncing back" of the spectrum from infinity after t = t* . Our considerations 
are inspired by Q. Exactly aXt — t* the spectrum becomes a pure power law, w"^'', at large frequencies. This can be 
understood as follows: any finite value of w, however large, is in the wake part of the self-similar solution when t — t* . 
This wake is a pure power spectrum. Therefore at i = i*, we have a well defined initial condition for the evolution 
equation. 

This spectrum is not a stationary solution (neither K-Z nor equilibrium) of the evolution equation. The subsequent 
evolution should follow the same principles as just before t = t*: the large frequency part of the spectrum has a 
typical timescale which goes to zero as uj 00. This timescale is the timescale for relaxation to a stationary K-Z 
spectrum with constant energy flux. Although the amplitude of this spectrum itself changes in the course of time, it 
does so more and more slowly after t = t* so that the changes induced in the K-Z spectrum by the change of energy 
flux become adiabatic relative to the infinitely short timescale for the large frequency evolution. This justifies our 
consideration of the large time spectrum as a K-Z spectrum although it is not strictly speaking stationary in time. 

Let assume now that the bouncing back of the K-Z spectrum from infinity is described by the same self-similar 
equation as before t = t* . One boundary condition is now different. The support of the self-similar spectrum now 
goes from zero to infinity and as a; — > 00, ~ uj~^° with xq being the K-Z exponent. Near small frequencies on 
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this stretched scale, the spectrum keeps the same power behaviom as before t = t* since the timescale there is long 
compared to the timescale t* — t &t the high frequency end. We expect that arguments which we shall present in 
sectionElare relevant here too, with the corresponding solution of the self-similar equation being unique and the free 
parameter required to adjust the trajectory reaching the low frequency behaviour along the stable manifold being the 
amplitude of the K-Z spectrum at infinity. The solution of the self-similar equation should then specify the shape 
of the "bend" in the spectrum seen in figure 01 and its time evolution will be then given be the same scaling of the 
frequencies in the similarity representation of the dynamical equations as we used before t = t* . Some further work 
will be required to make these statements more concrete. 



V. DERIVATION OF THE ANOMALOUS SPECTRUM FROM THE CRITICAL FRONT SPEED 

HYPOTHESIS 

In this section we present a heuristic derivation of the formula © promised in the introduction. The program is as 
follows. We first calculate the energy balance in the neighbourhood of the front tip. We will find that the condition 
that energy is conserved by the motion of the tip is equivalent to the condition, l|45|l . that the solution be self similar. 
Next we calculate the amount of energy, E{ti,t2), which enters the the frequency interval [uj* {ti),uj* {t2)] immediately 
behind the tip in the time interval [ti,i2]- The critical front speed hypothesis is the following : the physical system 
selects the unique value of the front speed such that limt^^f E{ti,t2) is finite and nonzero for all ti < t* . 

The basic energy balance equation for the front tip is as follows 

This can be expanded to read 



- + -]d.dt = , (58) 



'^'\E{iuM)~E{LoM)) doj^- [\piu;*it2),t)^P{ij*{ti),t)) dt (59) 

We observe that E{u!,ti) = since uj*{ti) < lo and P(ti>*(i2), t) — since Lo*{t) < L0*{t2). Hence we obtain the 
fundamental balance equation for the front tip, 

'^*(t2) rt2 

E{uj,t2)duj^ P{Lj*{ti),t) dt. (60) 
•(tl) Jti 

Let us now calculate P and E near the front tip using the expansion (|53() but keeping only the leading order term : 

F{r,)=ao{l-vf^^ + 0{{l-r,f/'). 

In terms of to and t, 

n{Lu,t) = aoc^*(i)"""^/^(^*(i) - w)^/' + 0{{Lu*{t) ~ w)^/^). 

We substitute this expression for n^{t) into equations and and keep only the leading power of uj*{t) — uj. 
For the energy, we obtain 

E{LJ,t) = ildUo LJ*(t)-"-3/2 w'i/" {LU*{t)-Ljf^^ 

+O{{L0*{t)-C0f/'). (61) 

For the flux we obtain 

8 

+0((u;*(i)-w)'/'). (62) 

In both these expressions we can write powers of lu as uj^ = (^*(^) ~ (^*(i) ^ '^))^ ^-i^d perform a Taylor expansion 
in {w*{t) — Lo). Since we are keeping only terms to leading order in {Lu*{t) — lo), we can consistently replace lo^ with 
uj*{t)y to yield the following expressions: 

E{LO,t) = VLdOQ LO*{ty^~^/^+'^^" {L0*{t)-L0f/^ 

+0{{uj*{t)-Lj)''/^). (63) 
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8 

+Oiiu;*{t)~iuf^^). (64) 

We now substitute these expressions into ()60|l and perform the integrals to leading order. Let us first consider the 
LHS: 

^'\{LU,t2)dLU = n^ao f '^^\*{t2y-'^'^+'^'°' {uo*{t2)~uof''^duj + h.o.t. 

= lndU0*{t2)-^-^'^+'"'' {U*{t2) - UJ*{tl)f'^ + h.o.t. (65) 

5 



Now consider the RHS: 



P{uo*{ti),t) dt - ^ flo / tj*(i)3(^o-^)-3/2 {uo* {t) - UJ* {ti)f '"^ dt + h.0.t. 



ti 



Let us change integration variables t — s- uj*(t): 

o "'w*(tl) V / 

We now use the self-similarity condition, (|45|l . to express the derivative in terms of w*(i) and perform the Taylor 
expansion trick again: for some power, y, we write : Lo*{t)y = {Lo*{tx) + (i^*(t) — ti'*(ti)))^ = uj*{ti)y + h.o.t. This 
gives 

= aijc^*(ii)3(-°-^)-3/2a;*(ii)-(3-o-2--d/a) / (^* - c^)3/2 do,* + h.o.t. 



(66) 



Equating H65() and H66|l and letting ti —>■ t2, we obtain energy balance criterion 



dt 

which is equivalent to the self-similarity condition that H45|l be time independent. 

Now suppose we take t2 — ti to be small but finite and allow t2 ^ t* . It is clear that the energy flux, (|65|l , entering 
the region [lu* (ti) , to* {12)] in this last increment of time before the front reaches infinity is either infinite, a finite 
quantity, or zero. Our hypothesis is that the flux is finite. It certainly cannot be infinite since the entire system has 
finite energy. It seems unreasonable that it should be zero since the front presumably requires a supply of energy to 
continue moving. A finite value for the integrated fiux requires that the power of w(i2) in H65|l is zero. This hypothesis 
leads to a unique value for the slope, which we shall denote by Xc- From l|65|) : 

lfldu;*it2r^-^/^+''^'' {Uj*{t2) - UJ*{tl)f'^ + h.o.t. 
5 

= ^lauj*{t2)-^-^'^+'"'^ {j^^ ' \{t2 h.o.t. (68) 

= c^*(i2)"^("-"°-^) ^J-{t2 h.o.t.. (69) 

5 

In order that this remain finite but non-zero as uj*{t2) ^ 00, we require 

X ^ Xc ^ xo -\ — (70) 

lia 

The physical intuition behind this hypothesis is the following. If x > Xc, the front speed is slower than the critical 
value. The power in (j69|l is positive and the energy in the tip diverges as w*(t2) 00. The front is moving too slowly 
for the amount of fiux fiowing into it so that energy begins to pile up at the tip. On the other hand, if x < Xc, the 
front moves faster than the critical value. Then energy in the tip decays to zero as w*(t2) 00 which means that the 
front is moving too fast for the amount of energy supplied to it. Physically we expect that the former situation would 
tend to speed up the front and the latter would tend to slow it down thus providing the system with a self-regulatory 
mechanism which selects the marginal slope, Xc- 
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VI. CONCLUDING REMARKS 

From the work presented here, we are beginning to get a clearer understanding of how the K-Z spectrum is set 
up in this model and the origins of the anomalous spectrum in the wake. The most obvious and important question 
which we have not addressed at all here is that of whether any of this analysis is relevant for the full kinetic equation, 
. It seems unlikely that the detailed structure of the nonlinear front would carry over to the integral version of the 
kinetic equation. Indeed, it is difficult to see how a sharp front tip could co-exist with the k-space integrations on the 
RHS of ((TJ. 

Nevertheless, anomalous behaviour has been observed ^ in numerical simulations of the full three-wave kinetic 
equation as already mentioned. This gives us reason to hope that some of the qualitative ideas contained here can 
be extended to more general kinetic equations. In particular, if the dynamics of the front is indeed regulated by a 
critical speed hypothesis of the type proposed in section then there is hope that an analogue can be formulated 
even in the absence of a sharp tip. 

One might also speculate that this mechanism might occur in the evolution of the Kolmogorov spectrum of hy- 
drodynamic turbulence since the direct cascade in this case is also of finite energy capacity. Unfortunately, in the 
absence of a closed kinetic equation for hydrodynamic turbulence it is not clear how one would even begin to address 
this question from a mathematical point of view. Nevertheless, for the purpose of stimulating debate, let us make 
the following conjecture. In far from equilibrium situations with stationary states which have finite capacity, the 
evolution towards the stationary spectrum takes place in two stages. In the first stage, < t < t* , which is rapid, the 
system attempts to close the connection with the dissipative sink at very high (or very low) wave- numbers. In that 
stage, entropy production is positive as the system attempts to explore all the available phase space subject to the 
constraint of energy conservation. The wake spectrum is steeper (shallower if the sink is at k = 0) than that of the 
final stationary state. This slope is determined by the requirement that a finite amount of energy per unit time is 
delivered to the front tip at all times less that t* . After t = t* , energy is no longer conserved but entropy production 
is still positive as the system now explores a larger volume of phase space. However, the entropy production now 
decreases as a new front with the K-Z spectrum in its wake (between the front and k = oo) moves towards lower 
wave-numbers and invades the steeper spectrum set up during the first stage of evolution. While it may be difficult 
to confirm these conjectures in a quantitative manner for the variety of situations to which we suggest these ideas 
apply, it should not be too difficult to establish them (or prove them incorrect) qualitatively. 
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Note Added in Proof 

We have subsequently analysed a second order model equation whose structure is similar to that of the differential 
kinetic equation studied in this paper but is analytically and numerically more tractable. We have found qualitatively 
similar behaviour. However the value of the anomalous exponent appears to differ from that which would be predicted 
by our critical front speed hypothesis. 
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APPENDIX A: ANALYSIS OF THE SELF-SIMILARITY EQUATION AND THE CRITICAL FRONT 

SPEED 

In this appendix we shall analyse the similarity equation H44|) further. We wish to check that it does indeed admit 
solutions which reproduce the critical behaviour which we have ascribed to the solutions of its antecedent PDE in 
section IVI Specifically, we are interested in solutions which match the F{ri) (1 — rj)"^^^ singularity at = 1 to a 
power law solution, F(r]) ~ ry"^ as 77 ^ 0. 

Let us write in the form 

(-^ + ) - <^(- - - -) i {"-^'i (?)) ■ 

and look for solutions of the form F{ri) ~ Arj^^ . Substituting this form we obtain 

i-y + x) ij^-'-y = 3A^y{y-l){xo-y)i3x„-3y-l) 

i2{x - xo) - t^d) 773-0-3^-2 (A2) 

We see that we can have the following solutions y = x where x takes one of the following values, 

X = (A3) 
X = 1 (A4) 
X — xq (A5) 
X = Xo- 1/3 (A6) 

(A7) 

The first pair are the thermodynamic spectra, the second pair are the K-Z energy and particle spectra respectively. 
A fifth special value of y is 

X = xo + Kd/2. (A8) 

What we observe from numerical solution of the O.D.E. 1A2|I . described below, is as follows. If we choose a value 
of X which is not Xc — xq + Kd/4:, then the solution, rj^^, makes a transition near 77 = to the state 77"^, where 
y — Xq + Kd/2. This behaviour near 77 = balances the leading order divergences on both sides of equation (|A2|) as 
77 — > 0. It might be noted, although it may have no relevance to the problem under consideration, that this value 
for the exponent, y, leads to front dynamics, (^*{t) — {t* — t)^ with zero b. This spectrum was never observed in the 
solutions of the P.D.E. 

To find solutions of (|44|) which do not exhibit the xq + Kdl2 scaling as 7; — > 0, we decided to perform a set of 
numerical experiments. Let us write out the RHS explicitly so that we can see exactly the equation which we wish to 
solve: 

' xF + ri— = -ifF^-— - 2sT]'-^F^ 



2{x -Xo)-Kd ' V dT] J ' drj'i ' drj^ 
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'dF\ fF 
drj J drf^ 

dF d^F 



/d^F 

+2s{s-1W-'f(^ 
\dr] 



(A9) 



drj drf 



The problem with integrating l)A9|l on a computer is that for generic initial conditions, the strong power law depen- 
dences of the right hand side on the independent variable render the numerics very susceptible to round-off error and 
numerical instability. For example, to study the system with 7 — 7/2, a = 1/2, d = 1, to which we gave a lot of 
consideration in section UTTl due to its relatively large anomaly, we are required to take 3xo + 2 = 22. 

To get around this difficulty, and to aid visualisation of the global properties of the equation, we make the following 
change of variables : 



Fill) 
dF 
drj 
d^F 
dif' 
d^F 
drf 



(AlO) 



where r = log(?7). By choosing 



I ( d \ Kd 



(All) 



we can cancel all the power dependence from the system and eventually recast IjAQjl as the the following autonomous 
fourth order system: 

dr 
dg_ 
dr 
dh 

> dk 



dr 



9 - 


af 




h- 


(«- 


1)5 


k- 


(a- 


2)h 






1 



(A12) 



2{x - Xq) - Kd 



(.9 + xf) - (a + 2s - 3)fk + 2fh^ 



~8g^h + Asg^ + Asfgh + 2s{s - l)fg^ ~ s{s - l)fh. 



Only the region, / > makes physical sense since n^^, and hence F, cannot be negative. This system is singular on 
the hyperplane / = 0. This system is much easier to integrate numerically and has the added advantage that we can 
determine the presence of fixed points which are not obvious in the original differential equation. Let us determine 
these points. It is obvious that the RHS of (|A12p has a trivial zero at (/, g, h, k) = (0, 0, 0, 0) = O but this is clearly a 
singular point due to the factor of on the LHS. A second pair of nontrivial (and nonsingular) fixed points can be 
shown after quite a bit of algebra to exist at the points P± = (/o, a/o, a{a — l)/o, a(a — l)(a — 2)/o), where 



/o = ± (18a(a + l)(a + {s ~ 2)/3)(a + (s - 3)/3))" 



(A13) 



We are naturally only interested in the point P+ since P- lies in the negative / region. The factors (a -f (s — 2) /3) 
and (a -h (s — 3)/3) are interesting for the following reason. If we substitute back in s = 3xo + 2 and the value of a 
from ljAll|l . these two factors are simply —Kdl2 and —Ki/2 respectively. At the transition points between finite and 
infinite capacity the fixed point runs away to infinity which suggests that it has a central role to play in organising 
the critical solution in the finite capacity case. 

Let us now look for a numerical solution of this system which mirrors the solutions of the differential kinetic 
equation. All the simulations in this section were done using the transformed system, ()A12|I . for which a standard 
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out-of-the-box adaptive Runge-Kutta routine seemed to work fine. Suppose we want a solution for the wake of the 
form F{r]) ~ r]~^ as 77 — > 0. In the new variables this is equivalent to demanding 



fiv) 
giv) 



x{x + l)/(r7) 
—x{x + l)(a; - 



as ?7 — > or as r — !■ 00. If a; < —a = a;o + which is the case, then we observe that the wake is described by a 

trajectory for which {f,g,h,k) — > (0,0,0,0) as r ^ —00. The wake is the singular point, O, of the system, (|A12(l . 
As we approach rj —> 1 we must reproduce the front structure described in section Hvl Thus we require a trajectory 
for which (/, /i, fc) — s- as r — s- 0, where F = (0, 0, cx), 00). We require a trajectory which links these two singular 
points. 

In our numerical simulations we rescaled the variables as follows, 



9 
h 

k 



I 

fo 
9 
afo 



a(a - l)/o 
k 



a(a-l)(a-2)/o 



This maps the point P+ — !■ (1, 1, 1, 1) for ease of visualisation. The only difference is that since a < 0, the signs of 
g and k are swapped so that the wake part of the solution must now approach O from the direction (O"*" , O"*" , O"*" , 0+) 
and the front tip, F, is now at (0"*", 0"*", 00, —00). 

It turns out to be difficult to find a trajectory linking _F — > O. We performed a series of experiments integrating 
backwards from 1 — ?/ = e, with e << 1, towards 77 = 0. Because we cannot specify initial data exactly at the 
singular point, F, we were required to manually tune the initial conditions quite a bit in order to to reproduce the 
^iv) '^-^ (1 — 77)^^^ structure near the tip. The remaining adjustable parameter is the equation is x. We again choose 
to study the case 7 = 7/2, a = 1/2, d = 1. Two generic types of trajectory emerge as we vary x. We visualise these 
trajectories in 4 dimensional phase space as a pair of projections of the actual trajectory onto the (/, g) and [h, k) 
planes respectively, hence the apparent intersections. 
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FIG. 7: Trajectory in {f,g, h, k) space for x = 7.56275, which is less than the critical value. 

FigureElshows what happens when x is slightly less than Xc- The trajectory leaves the front tip and heads towards 
the fixed point, P, but deflects to the right and is quickly attracted onto the line corresponding to F{ri) ^ ^ 
Figure IHl shows the corresponding trajectory when x is slightly greater than Xc- This time the trajectory deflects to 



17 




FIG. 9: Trajectory in {f,g,h, k) space for x = 7.56875 ~ Xc- Since it is impossible to specify the data exactly, if we integrate 
for long enough in r this trajectory eventually deflects away from O. 



the left as it approaches the fixed point, P, and is attracted towards the singular point, O. Before it reaches there 
however, it is deflected and is rapidly attracted back onto the F{ri) ^ gQi^^^ion again. The transition point 

between these two trajectories is illustrated in figure El Further analysis is required to determine the exact nature of 
the critical trajectory. It is not clear whether it will eventually deflect away from O for sufficiently small -q. Of course, 
in practical terms there is a limit placed on the extent of the scaling region by the left boundary of the inertial range 
which cannot extend all the way to in a real experiment. 

The system is very sensitive to the value of x. The values of x for the trajectories shown differ only in the 
fourth decimal place. In figure the form of the function F{r]) is shown for the three cases discussed above, after 
transforming back from the {f,g,h,k), r variables. The integration was again done for the case 7 — 7/2, a — 1/2, 
d = 1. Despite the small difference in the values of x in the equation the difference in asymptotic slopes is 7.56 for 
the critical slope versus 8.5 for the other two. 
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FIG. 10; Plots of the F{rj) associated with the trajectories in figures |7| |H| and |^ after converting from the {f,g,h,k) and r 
variables back to F and rj variables. The scales are log-log. 



APPENDIX B: OUTLINE OF NUMERICAL METHODS 

In order to study how the K-Z spectrum is set up from some given initial condition we wrote some code to solve 
(I13|l numerically. In this appendix we briefly outline the approach used. Write equation 1)1 3|l in the form 

-Q^ = f [nuj] + 9 [nj\ + h [n^] + F [nj\ - v \nj^ , (Bl) 

/K] = -K>^^'r? (B2) 
= -1k-Js^'~^n^ (B3) 



where 



2scj"-i 12(^1 + 2n 



duj ) doj dui'^ j 

-<'-)-^(-(£)^ <-) 

A„ = —uj--^ (B5) 
a 

s = 3x0 + 2. (B6) 

We are now including forcing and damping terms, F [n^] and v [rii^] which can be chosen as appropriate. We performed 
the following implicit time discretisation, 

it + At) - (t) d*n^ {t + At) d^n^ [t + At) 

At = — aJ3 — + 9M)] — g^, — 

+h [n^ {t)]+F [n^ {t + At)] -iy[nJt + At)], (B7) 

with the aim of enhancing the stability of the higher order derivatives. This can be rearranged to yield a time stepping 
algorithm in the form 

n^it + At) = [n^(<)] B [n^{t)] , (B8) 
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where 



L [nj\ 




(B9) 
(BIO) 



The time evolution operator, L [n^^], and source, B [n^^], are approximated using centred difference representations for 
the derivatives and a standard linear solver used to perform the inversion at each time-step. Implementing consistent 
boundary conditions is a tricky task. For all of the simulations presented here, the equation was solved from a 
compactly supported initial condition either in a frequency interval sufhciently large that the solution does not reach 
the boundary within the time allotted or with the damping chosen sufficiently strong to prevent the solution from 
ever reaching the boundary. 

Most of the simulations presented here are of a freely decaying initial energy distribution so F[ntJ\ = 0. The damping 
was chosen to be zero over most of the interval but increasing strongly for uj < lol and uj > ujr to produce regions 
of strong dissipation at large and small scales. We typically used approximately 10 grid points per unit frequency 
interval in our discretisation. The simulations presented here used up to 20000 grid-points. The number of grid 
points is practically limited by the fact that one must recompute and invert the time evolution operator, (|B9p . at each 
time-step. Lacking any reliable stability criteria for the nonlinear discretisation scheme described above, the choice of 
time-step was essentially made by trial and error. Once a time-step was found for which the evolution seemed stable, 
we ran with it. We were aided in choosing this time-step by monitoring the two conservation laws associated with the 
total energy and total particle number. Typically, and unsurprisingly, the numerics become unstable most easily at 
the front tip. This instability is worse for steeper spectra when the front travels faster. The necessity of reducing the 
time-step to resolve this structure as the tip increases in speed also put practical limits on the time interval we could 
simulate for. Retrospectively, we would like to do the numerics using an adaptive grid which tracks the front. Such 
a method would probably yield great increases in efficiency and stability. Unfortunately at the outset, we did not 
really appreciate that the code would be required to tackle this problem. We settled for validating our computations 
by checking that our final results remained unchanged when the time-step was reduced by a factor of two. 



